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Abstract 



Employing the path integral approach, we calculate the semiclassical equilibrium 
density matrix of a particle moving in a nonlinear potential field for coordinates near 
^ ' the top of a potential barrier. As the temperature is decreased, near a critical tem- 

' perature Tc the harmonic approximation for the fluctuation path integral fails. This 

is due to a caustic arising at a bifurcation point of the classical paths. We provide 
fS| I a selfconsistent scheme to treat the large quantum fluctuations leading to a nonlin- 

ear fluctuation potential. The procedure differs from methods used near caustics of 
the real time propagator. The semiclassical density matrix is determined explicitly 
for the case of asymmetric barriers from high temperatures down to temperatures 
somewhat below Tr. 



D^. 1 Introduction 

> 

;_j ■ Processes hindered by potential barriers play a central role in many areas of 

^ physics and chemistry When studying these systems, as a first step one 

has to determine the static behavior, that is the equilibrium density matrix 
Pi3 oc exp{—f3H). In barrier penetration problems the barrier height of the sys- 
tem is often large compared to the quantum mechanical energy level spacing. 
Then, a semiclassical evaluation of the equilibrium density matrix is justified. 
A consistent scheme to determine the coordinate representation of p/j in the 
semiclassical approximation is provided by the path integral approach 
Formally, the equilibrium density matrix may be viewed as a propagator in 
imaginary time. Hence, one first has to investigate the classical mechanics in 
the inverted nonlinear potential —V{x). Afterwords, the quantum fluctuations 
about the classical paths are determined. Since the classical equation of mo- 
tion, following from Hamilton's principle of least action, has to be solved with 
fixed endpoints, one finds in general a set of classical paths contributing to 
the path integral. If the classical paths, i.e. the minima of the classical action. 
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are well separated in function space, the contribution of the quantum fluctu- 
ations around each path is given by the simple semiclassical approximation 
where the fluctuation path integral is evaluated in the Gaussian approxima- 
tion. However, specific divergences, known as the problem of caustics, arise 
near those points where new classical paths become available. For a system 
with a potential barrier one encounters a caustic when the temperature is low- 
ered, specifically in the temperature range where thermally activated barrier 
crossing changes over to quantum tunneling ||^. In this region one has to go 
beyond the Gaussian approximation for the fluctuation path integral. For the 
case of a symmetric potential field, the semiclassical approximation near the 
caustic was investigated in a previous paper ^j. Here, we extend the results 
to asymmetric barriers. 

The paper is organized as follows. In section 2 we give a brief introduction 
to the path integral representation of the density matrix and the semiclas- 
sical approximation. In section 3 we then determine the classical paths and 
the classical action of a particle moving in the inverted barrier potential with 
endpoints in the vicinity of the barrier top. In section 4 we study the quantum 
fluctuations and determine the semiclassical density matrix at high tempera- 
tures and for temperatures near the caustic. Finally, in section 5 we give some 
explicit results for a barrier potential with cubic and quartic anharmonicities 
and present our conclusions. 



2 Path Integral Representation of the Density Matrix 
and Semiclassical Expansion 



The coordinate representation of the unnormalized equilibrium density matrix 
of a particle of mass m moving in a potential V{x) may be written as 



ppix,x') = Jv[x] 



^SeIx] 
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where the functional integral is over all paths x(r), < r < hj3 with a;(0) 
x{h(3) = x' . Each path is weighted by its Euclidian action 



up 
dr 







-mi;^ + V{x) 



(2) 



Since the density matrix can be regarded as the analytic continuation of the 
real time propagator to imaginary times t —ih(5, the representation (|1]) is 
frequently called the imaginary time path integral. The Euclidian action (0) 
describes the real time motion of a particle in the inverted potential —V{x). 
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Evaluation of the path integral (|l]) in an asymptotic expansion for h 
provides a consistent scheme for a semiclassical approximation. Thereby one 
first determines the maximum of the weighting factor, that is the minimum 
of This is given by the classical action 5'e[xci], where Xd is the clas- 

sical path solving the classical equation of motion following from Hamilton's 
principle (5S'£;[x] = 0. An arbitrary path in (|l|) reads 

x{t) = Xci{t) +y{T), (3) 



where ?/(r) describes the quantum fluctuations about the classical path. The 
fluctuations have to fulfill the boundary conditions y{0) = y{h(3) = 0. Using 
d^) the full action is then expanded around its minimum according to 




with the Frechet derivatives 
(5"5'£;[xci, y] = Jdri--- drn 





Sx{ti) ■ ■ ■5x(r„) 



yin)---yiTn 



(5) 



This way the dominant term (for ^ — > 0) is separated off and one is left with 
a functional integral over closed paths. In the simple semiclassical approxi- 
mation, the sum is truncated after the quadratic term leading to an exactly 
solvable Gaussian path integral pj-Q. If there exists a set {x"i} of classical 
trajectories in V{x), the above procedure must be performed for each and 
all contributions are summed to yield the semiclassical density matrix 

pMx,x') = E^e-^''"''Sl, (6) 

a V '-'a 



where Jq, = det{6'^SE[x]/Sx{Ti)6x{T2)\x=x'^^} is the determinant describing the 
Gaussian integral over the quantum fluctuations Clearly, the result in 

(H) is exact if V{x) is harmonic. Ja is given by the product of the eigenvalues 
A" of the second order variational operator 6'^Se[x]/Sx{ti)6x{t2)\x=x"^ as 

j„ = ivM!^ n K, (7) 

m 



where N is an appropriate normalization constant. As long as the second or- 
der variational operator is positive definite, i.e. A„ > for all n, the Gaussian 
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approximation gives the leading order fluctuation term for h ^ 0. But a prob- 
lem arises if one of the eigenvalues A„ tends to zero, e.g. as the temperature 
is lowered. Then, the quantum fluctuations of this mode become arbitrarily 
large and the simple semiclassical approximation breaks down. Generally, the 
vanishing of an eigenvalue defines a point where new minimal action paths 
in the potential V{x) become possible. This is well-known as the problem 
of caustics. In those regions where a caustic arises one has to go beyond the 
simple semiclassical approximation and has to consider higher order terms in 
the expansion of the action. 

3 Classical Mechanics near the Barrier Top 

To determine the density matrix in the semiclassical approximation we first 
have to study the classical mechanics of the system in the inverted potential. 
Hence, we have to solve the classical equation of motion with the boundary 
conditions x(0) = x, x{hj3) = x', where the endpoints x and x' are in the 
barrier region. A general barrier potential may be written as 



where the a„ are dimensionless coefficients that are of order 1 or smaller and 
Xa is the characteristic distance from the barrier top at which anharmonic 
terms of the potential become relevant. We assume that for coordinates in 
the vicinity of the barrier top the potential (H) is only weakly affected by 
anharmonicities. This means that the length scale Xa is much larger than the 
quantum mechanical length scale 



which is the variance of the ground state of an harmonic oscillator in the 
inverted potential. Correspondingly, 



is a small dimensionless parameter which will be serve as an expansion pa- 
rameter for the semiclassical approximation in the following. This investigation 
becomes more transparent if we introduce the dimensionless coordinate 

q = x/xo, (11) 







(9) 
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and the dimensionless time and inverse temperature 



a = OUT, 

e = ujh(3. 



'12) 



In terms of these dimensionless variables the scaled action S = Se/^ reads 



da 



(13) 



To determine the classical paths in the time interval < a < 0, it is convenient 
to use the Fourier series expansion 



-j^ oo 

^(^) = ^HQk sin(z/fccr) 
^ k=i 



(14) 



with the frequencies 
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(15) 



The above series continues the path outside [0, 6] as an antisymmetric and 
periodic path with period 26. The continued path has jump singularities at 
the endpoints of the interval [0, 6]. This artefact of the Fourier series expansion 
must be taken into account when calculating time derivatives of q. Hence, in 



([T3| ) one has to insert q^eg instead of q where 

greg(cr) = q{cr) - Ci : 6{a, 26) : -Ca : 5{a + 6, 26) : . 



(16) 



Here, the coefficients ci, C2 are determined by the jumps at 2k6 and {2k — T 

as 



ci = g(0+)-g(0-) = 2g 
C2 = g(^+)-g(r) = -2g' 



(17) 



and 



fc=— oo 



is a periodically continued (5-function. Inserting ([M]) and (|T^ into (|T3D, the 
action is found to read in terms of the Fourier amplitudes 



5 



1 

46 



(q - q') 



l\2 



1 oo 

-Y.[^kQl-2h{q,q')Qk + 4 

k=l 



+ g ^ ("£ ' " 

n=3 fci...fc„=l 



(-i)V 



where we have introduced the coefficients 
1 



z. n 

2 / dx ]^ sin (vr/c 



■jX) , 



(19) 



(20) 



and 



(21) 



as well as 



h{q,q') 



2'Kk 



g-(-l)V 



f22) 



Requiring that the variation 6S[Qk] = 0, we get the Fourier representation of 
the classical equation of motion 



AfcQfc — bk — O-n 



oo / g X n-2 oo 



n=3 



^kik2---kn-ikQki ' ' ' Qkn-i- 

fcl,---,fen-l = l 



(23) 



In general, the nonlinear equations (^) cannot be solved exactly. However, for 
small e and not too large 6 one may solve them pertubatively in a selfconsistent 



manner. 



3.1 High Temperatures 



Since the boundary values q, q' are in the barrier region and assumed to be of 
order 1 or smaller, for high temperatures, ^ <^ 1, the hk{q, q') are of order 1/9 
or smaller. The anharmonic terms in (|23| ) are then at most of order eO and 
the classical equation of motion reduces to 



XkQk = bk + {ed) 



(24) 
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Neglecting the corrections, we obtain the Fourier representation of the minimal 
action paths of an harmonic oscillator. When is inserted into ([ToD , the 
corresponding action is found to read 

S{z, r) = tan (^-) + ^ cot (^-] + O {eO') , (25) 



where we have introduced dimensionless sum and difference coordinates 

r = — - — , z = q — q . {2b) 



This approximation is only valid for endpoints within the barrier region and 
temperatures where the amplitude of the classical path remains within this 
domain. Clearly, the harmonic approximation fails if one of the coefficients 
Afc vanishes. When the temperature is lowered Ai reaches zero at the critical 
temperature 

= ^. (27) 



i.e. 6c = TT, and the corresponding amplitude Qi diverges. Corresponding 
singularities of the harmonic approximation arise for all temperatures where 
one of the frequencies Uk = 1. This is not surprising, rather it reflects the 
typical behavior of a harmonic oscillator in real time. For a harmonic potential 
the oscillation period 27r is independent of the amplitude of the path. Hence, 
for times 6 = nvr, corresponding to multiples of half the oscillation period, 
there are infinitely many paths with arbitrary large amplitudes connecting q 
and q' = (— 



3. 2 Classical Paths Near 



When the temperature approaches Tc from above, the eigenvalue Ai tends to 
zero and the harmonic approximation fails. Then, the anharmonic terms in 
(|23| ) become important. In general, the equation of motion cannot be solved 
analytically near T^. However, the classical paths and their actions can always 
be computed numerically. Here, we consider suitable potentials (see below) 
where only the amplitude Qi increases. For the bk{q, q') are now of order one 
or smaller and for perturbatively small e we shall assume that for tempera- 
tures near Tc the amplitude Qi may become at most of order e~^/^, while all 
other amplitudes are of order 1 or smaller. This order of magnitude of the 
Fourier coefficients will be confirmed below for appropriate barrier shape. The 
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A, 



Fig. 1. The Ai-r-plane is divided by the two curves r±{Xi) (sohd hnes) into two 
regions in which the cubic equation ( ^0[ ) has one or three (shaded region) solutions. 

equation of motion (|23|) then reduces to 

XkQk = bk- a^^-DukQl - a4 (0 ' DuikQl + O (e^/^) . (28) 

Within this approximation the amplitudes Q„ for n > 1 are coupled to Qi 
only. We have to solve a cubic equation for the mode amplitude Qi and linear 
equations for the other mode amplitudes Qn- To derive ( pSf ) we have made 
two assumptions. First, the parameter 04 is assumed to be positive and of 
order 1. In fact, the lenghtscale Xa in (|) can always be chosen such that 
a4 = 1 unless 04 is not positive. Second, the coefficient 03 can be taken positive 
without loss of generality. The coefficients a2„+i, n > 1 are assumed to be 
smaller than e^^^ which means that the barrier potential should be only weakly 
asymmetric. Only then, the cubic equation (|28| ) has real solutions Qi at most 
of order e~^/^ for all endpoints in the barrier region as assumed in deriving 
(p8| ) and correction terms of order e^^^. For larger asymmetries global features 
of the potential become relevant already in the vicinity of T^. To make the 
e-dependence more transparent we set 

Q = ^n^^'Ql- (29) 



Then, using Dm = S/Svr, -Dim = 3/4, and bi = Anr/O, the cubic equation 
for the mode amplitude Q takes the form 

3a^Q3 ^ ^ase-'/'Q' + \,e-'''Q = ^ r. (30) 
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Fig. 2. Bifurcation scenario for the classical paths for various boundary values r, a) 



r < Tc, b) Tc < r < 0, c) r = 
stable (unstable) solutions of 
and e = 0.01. 



0, d) r > 0. The solid (dashed) lines represent the 
The potential parameter are as = 1/5, 04 = 1 



From a mathematical point of view ( pO]) describes a pure (r = 0) or a perturbed 
(r 7^ 0) one-dimensional bifurcation problem with bifurcation parameter Ai 
M. The cubic equation ([30|) implies two curves 



r±(Ai 



3^-2±2(l 
Af- \ Af, 



3/2' 



(31) 



in the Ai-r-plane where bifurcations occur (Fig. p. Here, the critical point 
(Ac, is given by 



8046' 



2 



^3 



" 2n 
where 



Qi (32) 
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is the corresponding amplitude. Note that Ac = (4/3)''a|/a47r^ is independent 
of e. Further, Tc < 0, and Qc is of order one or smaller for 03, 04 > as assumed 
above. 

For Ai > Ac there is only one solution Q of (|30D which has the same sign as 
r. As the temperature is lowered, i.e. Ai is decreased, several cases must be 
distinguished. First, for r < Vc the high temperature solution continuously 
extents to the region Ai < Ac and the amplitude \Q\ grows (Fig. ga). Two new 
solutions with positive sign emerge when the temperature reaches the point 
where r_(Ai) = r (Fig. |T|). One of the new solutions is stable and the other 
unstable. The amplitude of the stable (unstable) path increases (decreases) 
with decreasing Ai. Second, for Tc < r < one first reaches a point where 
r+(Ai) = r when the temperature is decreased (Fig. |I|). At this temperature 
two new branches with negative sign but larger amplitude \Q\ appear (Fig. 0b). 
Again one path is stable and the other is unstable . With decreasing Ai the 
amplitude |Q| of the unstable path decreases and approaches the stable solu- 
tion which extends from the high temperature region. These two paths vanish 
when the temperature reaches the value where r_(Ai) = r. For temperatures 
below this point there exist only one stable solution until we arrive at the 
second solution of r_(Ai) = r. Here a stable and unstable path with positive 
sign emerge where the amplitude of the stable (unstable) path increases (de- 
creases) as the temperature is lowered further. For r = (Fig. |^c) we have a 
double point bifurcation scenario described below. Finally, for r > (Fig. |d) 
two new paths appear at the point where r+(Ai) = r. Again one solution is 
stable and the other one is unstable. 

In particular, for r = the solutions of (^Dj) are (Fig. ^) 



with Ac and Qc from (|32| ) and (^), respectively. One can show, and it will 
emerge below, that the trivial solution is stable for Ai > 0, the second solution 
is stable for Ai < Aq, while the third solution is stable for Ai < 0. 




(34) 



where Qo and Aq 



are given by 




(35) 
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3.3 Classical Action Near 



Having determined the classical path we are able to calculate the minimal 
action. Since (0) may easily be solved numerically, the action will be given 
as a function of the amplitude Q of the large amplitude mode. To get explicit 
values for the action the numerical value for Q must be inserted. Now, inserting 
( [I^ ) with the Fourier amplitudes determined from ( pS] ) and (|30| ) into (|T9|) and 
expanding in powers of e, we obtain after some algebra 



_ 4a3^(2 + 3A)r ^ _ 127^47- 3 _ 2aleV ^ 
3(47r2-^2),i/3^ 20A3 ^ e2/3 ^ 



^3 ^^3 

Here we have introduced the coefficients 



A = — tan - 



2 \ 2 / 7r2 - ^ 



2 



i Z = - cot - 
8 I2, 



4 47r2 - ^2 + 1^4^2 _ ^2 j ^^^1^2^ 97r2(7r2-e2)' 
which remain finite in the limit 6 ^ tt, i.e. 



r A 3 

hm A = - 

e^TT 2 

lim = 

limr = — (38) 

The first three terms in (|36|) reduce to the usual harmonic action for tempera- 
tures well above T^. Indeed, for Ai larger than order e^/^ the anharmonic terms 



11 



in (|30|) can be neglected and we have 



and hence 



-r^ tan 



O 



g4/3^ 



(40) 



While the harmonic action diverges for 6 ^ it, i.e. Ai — > 0, the full action (|36|) 
remains finite due to the contributions of the anharmonic terms. 

In the special case q = q' = 0, i.e. r = x = 0, one sees that for high tempera- 
tures the minimal action vanishes. With decreasing Ai two new paths emerge 
at Ai = Aq. At the bifurcation point the amplitudes of the new paths coincide, 
and to leading order the action is given by 



SciiQo, 0, 0)a,=Ao = ^a4g^e-2/3 + O (1) . (41) 

Since 5'ci(Qo, 0, 0)ai=Ao > 0, the trivial solution {Q = branch) is absolutely 
stable. At lower temperatures, i.e. smaller Ai, the solutions separate. Denoting 
the stable and unstable branch by Qs and Qu, respectively, the new amplitudes 
reach Qg = 2Qo and = 0, at Ai = 0, and the corresponding actions are 
given by 



S,i {Qs, 0, 0)^^^o = -49a^Qy^/^ + O (1) 

5ci(Q.,0,0),^^o = 0. (42) 
Now, ^ci(a,0,0) 

Ai=Ao < 0) s-i^d the Qg branch is absolutely stable. Between 
Ai = Ao and Ai = at the point Ai = 2Ac/3 the two stable paths have to 
leading order the same action and they exchange global stability. 

This example demonstrates the change of stability of the classical paths. At 
high temperatures there is only one minimum of the action in function space. 
With decreasing temperature a saddlepoint emerges at Ai = Aq with an action 
larger than the minimal action. This saddlepoint splits into a maximum and a 
second minimum of the action. The original minimum and the new minimum 
exchange global stability at Ai = 2Ac/3. Finally, the original minimum and 
the maximum meet in a double point at Ai = 0. 



For finite q and q' the cubic equation (|30D may easily be solved numerically. 
Again new extrema of the action functional arise as the temperature is lowered. 
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However, a change of global stability only occurs for small values of r. A more 
detailed discussion follows in section 4.4. 



4 Quantum Fluctuations 

Based on the classical paths and the associated extremal actions, we proceed 
to determine the equilibrium density matrix by expanding the action about 
the classical paths according to (|). 

We split an arbitrary path q{a) into the classical path qci{cr) and a fluctuation 
y{cr). The relevant fluctuations give a contribution to the full action of order h. 
Therefore, the semi classical expansion is only consistent if the classical action 
is also determined at least to order h that is to order 1 in the dimensionless 
units used above. In (^) we have neglected contributions to the classical action 
that are smaller than 1 for coordinates of order 1. On the basis of this result we 
may determine the semiclassical density matrix in the vicinity of the barrier 
top. For coordinates that are larger than order 1 a semiclassical expansion 
is of course still feasible, but mostly one has to use numerical methods to go 
beyond the approximate results derived in the previous section. The exception 
are some special potentials, for instance the Eckart barrier potential, where 
the classical mechanics can be solved exactly 

4-1 Expansion about the Classical Path 

To evaluate the pathintegral an arbitrary path is decomposed into 



and one has to expand the action about qd. If there exist a set gcil*^) of stable 
classical paths, the action is expanded about each Then the density 

matrix (0) may be written as 



where Sc\[Q, z, r] is the scaled minimal action (|T3D and Z denotes an appropri- 
ate normalization constant. f{Qa, z, r) is a functional integral over the paths 
y{a) given by 



(43) 



Pe{z,r) = - Y,(iW{-Sci[Qa,z,r]) f{Qa,z,r), 



(44) 




(45) 
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where S"^S[qci,y] are scaled Frechet derivatives according to @. Since the 
fluctuations t/(cr) have to satisfy the boundary conditions y{0) = y{6) = 0, we 
use the Fourier representation 

-j^ oo 

y(^) = nJl^k sin(z/fccr). (46) 

fc=i 



Then, the integration measure becomes 



Jm- -nlWs^/civ;-), (47) 



-oo 



where can be determined from the limit of a free particle, i.e. uj ^ 0. With 
the Fourier representations of the classical path (0) and the fluctuations (^61), 
respectively, one finds for the scaled Frechet derivatives 

1 oo 1 oo ^ s n-2 

^^%c,..] = ^gA,n^ + ^Ea„(n-l)(^) 

OO 

X J2 Dk,...k,,Qk,---Qk„^,Yk^^,Yk,^ (48) 

and 

xmcr 1 1 ^ {n-l)\ /e\"-2 

oo 

X X! ^ki...k„ Qki ■ ■ - Qk (49) 

for m > 3. Here, the coefficients and -Dfci...fc„ are defined in ( pTD and (pQ]), 
respectively. The first term of (^Sf ) arises from the bilinear terms of the action, 
while the remaining terms are due to the anharmonicity of the potential. 
Clearly, from ( ^9]) we see that the semiclassical expansion proceeds in powers 
of e. Of course, the fluctuation integral ( ^5| ) cannot be done exactly. However, 
for small e we may evaluate the functional integral (P5| ) perturbatively. To 
do so, we first have to solve (|23|) perturbatively for small e. In particular, 
this yields the order of magnitude of the Fourier coefficients Qk that depends 
on the temperature range. Afterwards, one estimates the size of the Fourier 
coefficients of the relevant fluctuations. The exponent of the integrand of 
(|45|), given by (^Sf ) and (|49|), is an expansion in powers of Y^ with coefficients 
depending on the parameter e, if we express the order of magnitude of the Qk in 



14 



terms of e. In a first step the size of the is estimated from the assumption 
that the terms in the action quadratic in 1^ are of order 1. However, this 
estimate is only correct if the terms of higher order in vanish in the limit 
e — >■ 0. If there are higher order terms that grow in this limit one gets a new 
(smaller) estimate of the Yk by requiring that the leading order stabilizing 
terms are of order 1. Again the estimate depends on the range of temperatures 
considered. 



4-2 High Temperatures 



Above, when studying the classical mechanics we have seen that for tempera- 
tures 9 <^TX and for coordinates near the barrier top the anharmonic terms are 
neglectible. In the limit of small e and for high temperatures the magnitude 
of the Fourier coefficients Qk is of order 9. Assuming that the coefficients 
are of order 6'^/^ or smaller we see that the first term on the right hand side 
of (^) is of order 1, while the remaining terms of (|48|) and (|^) vanish in the 
limit e — > 0. Therefore, for the relevant fiuctuations the exponent of (^) is 
given by 

-| oo 

S[q]-S[q.,] = -Y.^nY^ + 0{e9'). (50) 



Thus, the path integral 
integrals 



reduces to a product of independent Gaussian 



-| oo 



— / dYfeCxpf — ^AfcFi? 



(51) 



where the eigenvalues of the second variational operator are given by (^I]). 
Now, the fiuctuation prefactor for the equilibrium density matrix is obtained 

as 

/ = [47rsin(0)]"'/' (52) 



where the normalization in ( |5lD is derived from the free particle limit. Combin- 
ing the action (|25|) and the fiuctuation prefactor (|5^ , the equilibrium density 



matrix for coordinates near the barrier top and for high temperatures emerges 
as 



Pe{z,r) 



1 



1 



47rsin(^) 



exp 




(53) 
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Hence, in this temperature and coordinate range, we regain the well known 
result for the inverted harmonic oscillator. Clearly for 6' — > vr, this result 
becomes singular. Then, the mode Yi undergoes large fluctuations and one 
has to go beyond the simple semiclassical approximation. 



4-3 Quantum Fluctuations Near Tc 



Since the eigenvalue Ai of the second variational operator vanishes for T = Tc, 
the simple semiclassical approximation fails near Tc. Under the assumptions 
made in section 3.2 about the potential parameters a„ we were able to solve the 
classical equation of motion ( PBD near the critical temperature perturbatively 
for small e. In this temperature range all amplitudes Qk for k > 1 are of order 
1 or smaller, while the magnitude of the amplitude Qi is of order e~^/^ or 
smaller. We now estimate the magnitude of the Fourier coefficients Yk. The 
amplitudes Y^. for > 1 are of order 1 or smaller since all eigenvalues for 
k > 1 are of order 1 in this temperature range. Hence, the fluctuations Y^ for 
k > 1 are bounded by the first term on the right hand side of (^81) . 



Near the critical temperature the amplitude Yi can become much larger than 
1. To estimate the size of Yi we first consider the terms of (^) that are 
quadratic in Yi. With (|33) and (p3|) we have 



-| oo 
8^ k=2 



1 

8^ 



Ai - A, + 9a4e2/3 (g - Q^f) 



(54) 



From (153) we see that the coefficient 
Ai(g) = Ai - A, + 9a4e2/3 (g - Q^f 



(55) 



vanishes at the critical point, where Ai = Ac and Q = Qc- Thus the size of Yi 
is only bounded by higher order terms in Yi. If we assume that Yi is of order 
g-i/2 Qj- smaller near the critical point, the fourth order variational term gives 

S'S[qcuy] = ^,e'Y,^ + o{l). (56) 



Indeed, this term limits the size of the mode amplitude Yi near Ai = Ac to the 
assumed order of magnitude. With this estimate for Yi and in a narrow range 



16 



around the critical point, where Ai(Q) < (9(e) one gets for the action of the 
fluctuations 



S[q] - = - E >^nY^ + ViQ^ Yi) + 0{e 

n=2 



1/2N 



(57) 



with the fluctuation potential 



1 

8^ 



(Ai - Ae + 9046^/3 (Q - g,)') 



£4/3 



+3a4— (Q-Qe)n' + ^(^) 



3a4 /e^ ^ 



(5J 



From ( p5D we see that the coefficient Ai((5) of the harmonic term may van- 
ish. This is indeed the case for values of r and Ai on the two curves r^(Ai) 
introduced in (|3T|). When these curves are crossed, the classical equation of 
motion (|30|) has a bifurcation. Near the bifurcation the remaining terms of 
the fluctuation potential always constrain Yi to fluctuation amplitudes of or- 
der e~^/^ as assumed above. For Ki{Q) > the fluctuation potential has only 
one minimum at Yi = 0, corresponding to one solution of the classical equa- 
tion of motion. For Ai{Q) < the fluctuation potential has three extrema, at 
Yi = and at 



Y^ 



-36 

^2/3 



(Q - Qc) ± 



402 



3a4e" 



(Ai - \c) 



302 
7473 



(g - Qcf . 



(59) 



Specifically, let us discuss the fluctuation potential for fixed Q = Qc- Hence, 



according to (|30|), we consider values along a line r = re(Ai) in the Ai-r-plane, 
where 



re(Ai)=r, (3^-2). 



(60) 



For fixed Q = Qc the fluctuation potential has the simple form 



V{Qc,Y^] 



1 

86 



(Ai - Ae) Y,' + ^a, ) Y,' 



(61) 



which is symmetric about Yi = (Fig. |). According to (|59|) for Ai < Ac the 
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Fig. 3. The fluctuation potential (| 
extrema Y± are at 



/ 4^2 



(Ae-Ai). 



Y+ 
for Q = Qc and temperatures near Tc 



(62) 



For Ai — Ac > e the fluctuation potential has only one minimum Yi = and 
the second order variational operator constrains Yi to fluctuation amplitudes 
smaller than order e"^/^ this region the simple semiclassical approximation, 
i.e. the truncation after the quadratic term in the fluctuation potential, is valid. 
For Ai in the region < Ai — Ac < e, the curvature of the fluctuation potential 
at the minimum at Yi = becomes smaller and the fluctuation amplitude of 
the marginal mode is stabilized only by the quartic term. When Ai is decreased 
further, the minimum at Yi = becomes a maximum and new minima arise 
at Y± according to the bifurcation scenario discussed for the classical paths. 
Indeed, one can show that the extrema of the fluctuation potential V{Q,Yi) 
are the solutions of the cubic equation (^) for Q apart from the scaling factor 
26'e-2/3 introduced in (|2|). In the region — e < Ai — Ac < the new minima 
are not well separated by the local maximum at Yi = and fluctuations from 
one minimum to the other occur. Hence, in this region the fluctuations are 
also stabilized by the quartic term in the fluctuation potential. When the 
temperature is lowered further the barrier height of the local maximum of the 
fluctuation potential becomes larger and the fluctuations about the minima 
decrease. Then the system is stabilized near the stable classical paths and a 
Gaussian approximation for the integral over the fluctuations around these 
paths is again appropriate. 
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4-4 The Density Matrix Near Tc 



Well above the density matrix is given by ([53|). As is approached the 
fluctuations grow, and one has to take into account the full fluctuation poten- 
tial which stabilizes the fluctuations by the quartic term. In the vicinity 
of the critical point (Ac, Tc) the various solutions of the cubic equation (|30|) are 
not well separated and the question arise, which branch must be inserted into 
the action and the fluctuation potential. However, one can easily show that 
the path integral over the fluctuations is independent of the branch chosen. 
Expanding the action around the classical paths one is left with the fluctuation 

I) it may be written as 



path integral (|45|) . Using (]47|) and 



/W=n(^\/^/dl^.)exp(- 



-| oo 
n=2 



Ai 



Att sm{9) 



K{Q), 



(63) 



where 

K{Q) 



J dY^exp{-V{Q,Y,)) 



(64) 



gives the contribution of the marginal mode Yi. Now, with (|36|) and ( |64|) the 
equilibrium density matrix in the vicinity of the critical point (Aic, Tc) reads 

Peiz, r) = ^^/^^5^^(<5) exp {-S^{Q, z, r)) . (65) 

Usually, explicit results for the density matrix can be calculated only numer- 
ically. After solving the cubic equation of motion ( pOD for given Ai and r one 
has to insert Q into the action (|36|) and the fluctuation potential (|58|) and 
determine the fluctuation integral (0). All steps involve rather simple nu- 
merics only and some results will be presented in section 5. As we have seen 
in section 4.3 the full fluctuation potential is needed when the coefficient of 
the second order term Ai{Q) is of order e or smaller. For |Ai(Q)| larger than 
order e the Gaussian approximation for the fluctuation potential is appropri- 
ate. Moving away from the critical point (AcT^), where Ai(Q) = 0, we reach 
values of |Ai(Q) | of order e by varying Ai by AAi of order e or by varying Q by 
AQ of order e^/^. The corresponding variation Ar of the coordinate r is then 
determined by the cubic equation (|30D. In the flrst case, varying Ai for flxed 
Q = Qc, we get Ar = O {as). In the other case, varying Q for flxed Ai = Ac, 



19 



Fig. 4. Regions in Ai-r-plane in which various results for the density matrix are vahd. 
Around the critical point (Ac, rc) in the dark region result (65) has to be used. In the 
white region, where only one solution of (30) exists, result (67) is appropriate, while 
in the shaded region result (70) holds. Along the dashed line re(Ai) introduced in 
(60) one has Q = Qc 

the order of magnitude of the variation of r is Ar = O (e^^^) • Since Q = Qc 
along the line re(Ai) introduced in (|60|), the region around the critical point 
where the usual WKB approximation breaks down is of the form of the dark 
rectangle shown in Fig. 4. In this region the result ( |65|) must be used. 

For Ai{Q) > e the fluctuation potential can be simplified to read 

ViQ,Y,) = ^A,iQ)Y,' + oil), (66) 



where the terms omitted are smaller than order 1 since the fluctuations Yi are 
now smaller than e^^/^. K{Q) is therefore a Gaussian integral and the density 
matrix fB5|) reduces to 



For high temperatures, i.e. 6 <^ n, one sees from ( ^OD and (pSj) that the result 
([67| ) matches into the density matrix in the high temperature limit (|53|) . 



On the other hand, as discussed in section 3.1, for temperatures below the 
critical region the cubic equation has two stable solutions which we call here 
Qsi and Qs2 where Qsi < Qs2- As we have seen above, for Ai{Q) < —e these 
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two branches are well separated in function space, and a Gaussian approxima- 
tion for the fluctuations around these stable paths is appropriate. Choosing 
Q = Qsi the fluctuation potential has a minimum at Yi = and it can be 
expanded to give 

V{QsuY,) = La,(Q^,)Y,' + o{1), (68) 

while near the other minimum at Yi = Y^ — y_ the fluctuation potential takes 
the form 

V{Qsi, Yi) = SciiQs2, z, r) - SciiQsi, z, r) 

+ W^2) [Y, - (r+ - Y.)f + 0(1). (69) 
The density matrix therefore reads 



[z,r) = - 



Ai 



zv47rsin(^) \jx;(q:;~) 



exp{-SciiQ si, z,r)) 



Al{Qs2) 



exp{-Sci{Qs2,z,r)) 



(70) 



This result matches onto (^) for values of Ai and r where the coefficient 
Ai{Q) becomes of order e or smaller. When the difference of the actions of the 
two stable paths \Sci{Qsi, z, r) — Sci{Qs2, z, r)| is of order 1 or larger, only the 
path with smaller action contributes to the path integral. The density matrix 
then reduces to 

P'^'^'^ = W 4^sin(g^iA,(a) ^"P^-^"^^-"'^^^' ^''^ 



where Qs denotes the path with smaller action. 



5 Discussion and Conclusions 



We have studied the equilibrium density matrix of a quantum particle near 
the top of a general anharmonic potential barrier in the temperature range 
where large quantum fluctuations render the harmonic approximation of the 
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-2 2 

Fig. 5. The potential (8) with potential parameters as = 1/5, 04 = 1, and a„ = 
for n > 5. 

potential insufficient even for coordinates in the vicinity of the barrier top. We 
have shown how the density matrix can be calculated by means of the path 
integral technique employing an expansion about classical paths. In a narrow 
range around the critical temperature Tc and in the vicinity of the critical 
coordinate Tc the Gaussian approximation for the quantum fluctuations, i.e. 
the simple WKB expansion, fails and one has to take into account higher order 
variational terms. The anharmonicity of the potential remains essential below 
the critical temperature. 



To illustrate our results we have evaluated the diagonal part of the density 
matrix, i.e. the position distribution function P{q) = p{z = 0,r = g), for a 
system with the potential parameters 03 = 1/5, and 04 = 1 while all other 
coefficients a„ = 0. Thus the potential (|) reduces to 



_^ 2 X 1 / X 

15 Xa 2 V Xa 



(72) 



describing an asymmetric double well (Fig. |^). Using the dimensionless for- 
mulation (0) the dimensionless potential V = V/huj reads 

V{q) = -\q' + ^^q'+'^q\ (73) 



The explicit results given above hold for weakly asymmetric potentials with 
a| < e <^ 1. This is the case for e = 0.01. 
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Fig. 6. The position distribution Pgiq) at temperatures near Tc as a function of q 
near the barrier top. 

In Fig. ^ the equilibrium position distribution function P{q) is depicted for 
various temperatures near Tc. [The normahzation factor for all figures shown 
in this section is choosen as Z = 1.] For high temperatures the position dis- 
tribution function is given by the harmonic approximation that is symmetric 
around the barrier top at g = 0. Anharmonic terms in the equation of motion 
([30| ) become relevant with decreasing temperature. Accordingly, the minimum 
of the distribution is shifted and the distribution becomes more asymmetric 
as the temperature decreases. 

In Fig. 7 we show the diagonal part of the density matrix as a function of 
the eigenvalue Ai. At the barrier top, i.e. at g = (Fig. 0a) and for high 
temperatures there exist only one solution of (|30D , Q = 0, and the distribution 
function is determined by (0). At Ai = Aq a new stable solution Qs of (|30|) 
emerges. Between Ai = Aq and Ai = it exchanges global stability with the 
Q = solution and for lower temperatures the position distribution function 
is given by (|7lD . While the results (^) and (^) diverge at Ai = and Ai = 
Aq, respectively, the general formular (^) matches between these results and 
remains finite. 

For the critical value of the coordinate q = and at Ai = Ac there exists only 
one solution Q = Qc oi (^). The simple semiclassical approximation (|67|) 
diverges at Ai = Ac and we have to take the general result (^). In Fig. 7 we 
see the matching of the result (^Sj) with the result (|67[) vahd below and above 
the critical region. 

In summary, we have shown how to evaluate the semiclassical density matrix 
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Fig. 7. The position distribution Peiq) at temperatures near Tc as a function of the 
eigenvalue Ai for (a) g = and (b) q = Vc- The soHd hue represents the general 
result (65). Furthermore, various approximation are shown. In (a) the dashed line 
represents result (71) while the dotted line represents result (70). The dotted-dashed 
line represents result (67), which is for q = the same as the high temperature WKB 
approximation (53). In (b) the dashed lines represent result (67). 



consistently near the critical region. For 03 = one regains the results for 
the anharmonic symmetric barrier potential investigated previously [^]. In 
this case the critical point (Ac,Tc) = (0,0). A comparison of the semiclassical 
density matrix with exact results will be given elsewhere. 
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